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In order to incorporate spatial inhomogeneity due to nonmagnetic impurities, Anderson proposed a BCS- 
type theory in which single-particle states in such an inhomogeneous system are used. We examine Anderson's 
proposal, in comparison with the Bogoliubov-de Gennes equations, for the attractive Hubbard model on a system 
with surfaces and impurities. 



The procedure for examining surface and im- 
purity effects on a microscopic level is by now 
well established. One uses a mean field-like de- 
coupling, with potentials which are determined 
from self-consistency requirements. These poten- 
tials are then used in the effective Hamiltonian, 
which is numerically diagonalized. This process is 
continued until self-consistency is achieved. This 
is the essence of the Bogoliubov-de Gennes (BdG) 
j|| formalism. 

An earlier proposal was suggested by Ander- 
son [jlj, in which the single-particle problem is 
first diagonalized. Eigenvalues and eigenstates 
are obtained, with which one can formulate the 
BCS problem, but in a vector space associated 
with these eigenstates. In certain situations, the 
single-particle problem can be obtained analyti- 
cally (open boundaries, for example p|), or it can 
be obtained numerically with significantly less ef- 
fort than required by the full BdG process. In 
these instances it would be advantageous to uti- 
lize the Anderson prescription. In this paper we 
report on some test cases to evaluate the Ander- 
son prescription. 

The BdG equations are well documented 
In this work we utilize the attractive Hubbard 
model, with open boundaries, and with the pos- 
sibility for single site impurity potentials. The 
resulting equations are M : 



E n u n (£) 



A u >u n {t') + A e v n (e) 



(1) 



E n v n {i) 



where 



A U 'V n {l') + A* e u n {£) 



(2) 



Aw = -t ^Sj^ e -s+Se',e+s) +&w yVj-fi+eij .(3) 

The self-consistent potentials, Ve, and A^, are 
given by 

A i = \U\^2u n {£)v*(e)(l-2f n ) 



v e = -\u\J2 



\u n (£)\ 2 f n + \v n (£)\ 2 (l-f n ) 



(4) 



, (5) 



where \U\ is the strength of the attractive inter- 
action, the index n labels the eigenvalues (there 
are 2N of them), the index £ labels the sites 
(1 through N), and the composite eigenvector is 



given by ^ v " J , of total length 2N. The sums in 

Eqs. (QJ^) are over positive eigenvalues only. The 
other physical parameters are the single-particle 
hopping, t, the single site impurity potentials, q, 
and the chemical potential, fi. The /„ is the Fermi 
function, with argument [3E n , where (3 = jr^y, 
with T the temperature. The single site elec- 
tron density, nt, is given, through Eq. (|^), by 
V t = -\U\?£. 

These equations are iterated to convergence, 
with results to be presented below. 

The alternative Anderson formalism jij first 
solves for the eigenvalues and eigenstates of the 
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'non- interacting' problem, i.e., 

E° n w n (l) = J2 A °u^n(i% (6) 
I' 

where 

A °w = -t^^e'.e-s + Se'^+s^ -5w ^/i-e^J.(7) 

The N x N matrix equation is solved for 
its eigenvalues E® and eigenvectors w n . This 
amounts to determining the unitary matrix Uf n 
that gives a basis for the electron operators 

n 

which diagonalizes the single-particle Hamilto- 
nian. From this matrix we determine the trans- 
formed electron-electron interaction: 

Vnm^n'm' = —\U\ U^ n U^ m Ul n >Ul m ' , (9) 

which now mediates the (generally off-diagonal) 
electron-electron interaction. The gap and num- 
ber equations are derived in the usual way; they 
are in general complicated — the gap is a function 
of the quantum label n and the chemical potential 
is shifted by an n-dependent quantity. Once these 
are obtained, we can transform back to real space, 
and examine the gap function or the electron den- 
sity, for example, as a function of position. 

Figure 1 illustrates the gap parameter obtained 
by the BdG formalism as a function of position, 
for all densities, in the case of open boundary 
conditions (OBC), in one dimension. Results in 
higher dimension will be very similar jq]. Vari- 
ations in the gap are strongest near the bound- 
aries, as expected, and the Anderson prescription 
is reasonably accurate in reproducing the oscilla- 
tions (not shown). In Fig. 2 we show the gap 
as a function of position for the case of a sin- 
gle impurity (at site 16) with a repulsive poten- 
tial, with periodic boundary conditions (PBC). 
As expected, the gap is suppressed at this site, 
and once again, the Anderson prescription semi- 
quantitatively reproduces the BdG result. 
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BdG : |U| / t = 2 N = 32 OBC 

A (i) / t 




Figure 1. A(i) for all densities n e (OBC). 




Figure 2. A(i) for n e = 0.9 (e(16) = 1.0 1). 
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